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Abstract 

We study power corrections to a recently introduced family of event shapes, 
the class of angularities, within the formalism of dressed gluon exponentiation 
(DGE). We find that the universal scaling rule for the leading power corrections 
deduced from resummation also holds when taking renormalon enhancements into 
account. The scaling is due to boost invariance of eikonal dynamics in the two- 
jet limit, which we recover in the context of DGE. Furthermore, dressed gluon 
exponentiation provides an ansatz for non-leading power corrections that violate 
the scaling. These non-leading corrections are further suppressed by non-integer 
powers of the hard scale. 
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1 Introduction 



Event shapes [H El IS] are infrared-safe generalizations of jet cross sections, and they 
are among the best suited observables to test our understanding of QCD. Event shapes 
distributions, in particular, describe the behavior of colored radiation in the final state 
of hard collisions in all of phase space. Thus, although infrared safe, they are sensi- 
tive to emission at all scales, and provide a unique tool to probe the interface between 
perturbative and non-perturbative QCD. 

A consequence of the sensitivity of event shapes to long-distance effects is that pertur- 
bative calculations are far from straightforward even at high energies, where the strong 
coupling is sufficiently small. Specifically, in the narrow-jet limit radiation is dominated 
by secondary partons which are either soft or collinear to the primary quarks emitted at 
the hard scattering. As a consequence, fixed order computations receive large logarith- 
mic corrections, which need to be resummed to all orders to obtain reliable quantitative 
predictions H El El 13 13 • 

In addition to these large perturbative corrections, there are corrections due to con- 
finement that cannot be treated within the perturbative expansion, and are suppressed 
by powers of the hard scale. Their main effect is to widen the distribution of radiation in 
the final state, shifting the peak away from the narrow jet limit. Although long-distance 
effects are intrinsically nonperturbative, the study of resummed perturbation theory can 
give valuable information about the size of their contributions to the cross section, as we 
will illustrate below. 

In the following we will consider event shapes in e + e~ annihilation. Power corrections 
to these observables have been the focus of intense research during the past several years 
(for recent reviews see for example [HIE])- It turns out that for mean values of such event 
shapes the effects of power corrections can be efficiently summarized in terms of a single 
parameter, essentially a shift of the perturbative distribution by an amount proportional 
to an integer power of l/Q, where Q is the overall center-of-mass energy [TmiTT)rr21 ll3|. 
It was found that these additive corrections are to a certain degree universal. This 
universality property can be used phenomenologically for precise measurements of the 
strong coupling from event shapes |14j . 

The case of differential distributions, however, is more complicated. Distributions 
probe scales that are even smaller than those which dominate mean values, thus sub- 
leading power corrections also need to be taken into account. For a large class of such 
corrections this can be done by introducing a nonperturbative shape function ^j. Re- 
summed perturbation theory, which displays sensitivity to soft emission at power ac- 
curacy, imposes constraints on models for these shape functions. Phenomenologically, 
the models are quite successful in describing the data, again with only a small set of 
nonperturbative parameters. Universality, however, is generically lost at the level of 
subleading power corrections, although certain classes of event shapes turn out to have 
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closely related nonperturbative behavior jTEJ [T71 [TH] d] . 

A class of event shapes of this type depending on a real parameter a was introduced 
in [20] • These observables provide an interpolation between the thrust [T|, corresponding 
to a = 0, and the jet broadening a — 1. In this paper, we will refer to event shapes 
in this class as angularities. As shown in |2TH [HJ 1221 122] > these event shapes can be 
studied analytically, as functions of a, both at the perturbative level, performing an 
all-order resummation of Sudakov logarithms, and at the level of long-distance effects. 
Using factorization arguments, it was found that leading power corrections to all the 
observables in the class of angularities are described by the same shape function, up to 
an overall scaling factor j2Sj- This remarkable property is closely related to the boost 
invariance of the soft radiation emitted in the two-jet limit. 

We will study power corrections to the distributions of angularities by using renor- 
malon methods. Specifically, we will use Dressed Gluon Exponentiation (DGE) [21] to 
build a model of the shape function. We will recover the scaling property discovered in 
[23]) and we will study the pattern of violations of scaling due to the contributions of 
collinear radiation. 

Dressed gluon exponentiation combines renormalon calculus with Sudakov resum- 
mation. On the one hand, taking into account single (dressed) gluon emission in the 
large /3o-limit, it identifies the leading source of factorially divergent behavior, character- 
istic of the asymptotic nature of the perturbative expansion. This is then upgraded to 
multi-gluon emission in the Sudakov limit by exponentiation of the single gluon result, 
employing as usual a Laplace transform to enforce momentum conservation. In the con- 
ventional resummed expressions, the asymptotic behavior manifests itself as a singularity 
in the perturbative running coupling at small scales, which introduces an ambiguity in 
the form of power corrections. The DGE formalism provides a definite prescription to 
deal with these ambiguities, and thus strongly constrains the form of nonperturbative 
shape functions. Dressed gluon exponentiation has been previously applied to several 
other event shapes, including the thrust and the C-parameter [TTl [THl HHI I^HJ , as well as to 
processes involving heavy quarks [26 , and inclusive hadronic cross sections such as DIS 
and Drell-Yan [231 12 7j. Here we will apply it to study power corrections to angularities. 

An interesting aspect of the study of angularities with this method is the fact that 
the scaling rule found in Ref. [23] was shown to be closely related to the boost invariance 
properties of the eikonal cross section describing soft emission in the two-jet limit. This 
property is not at all apparent in DGE, where in practice the single gluon cross section 
is computed with a "massive" gluon, along the lines of the dispersive approach. Since 
the introduction of the gluon virtuality breaks boost invariance, and since in addition is 
necessary to account for sub leading logarithms in the large rif limit, it is not a priori clear 
whether the scaling will survive. Here we will show how the effects of boost invariance are 
recovered in the Sudakov region. In the two-jet limit, only the logarithmically divergent 
bremsstrahlung spectrum contributes to the dressed gluon cross section, and the gluon 
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mass acts as an effective cutoff, with precisely the weight required to reconstruct an 
unweighted rapidity integral. 

We start by briefly reviewing the definition of the class of angularities, and the scal- 
ing rule for nonperturbative corrections found by analyzing the resummed expression. 
In Sees. El and 0] we will then construct a model of the shape function for angularities by 
means of DGE. Concentrating on soft radiation, we will recover the scaling of leading 
power corrections, and study how it arises in the context of a massive gluon calcula- 
tion. Finally, we will observe that DGE suggests an intricate, non-universal pattern of 
subleading power corrections arising from collinear radiation. All such corrections are 
however suppressed by non-integer powers of the hard scale, with a degree of suppression 
growing as the parameter a becomes large and negative, where the event shape becomes 
correspondingly more inclusive. Technical details are given in the Appendix. 

Although detailed phenomenological studies will have to be deferred to future work, 
when and if experimental data become available, we emphasize that such studies would 
be of considerable interest, and should be quite easy to perform for experimental col- 
laborations. The scaling rule found in j2H| and recovered here is in fact a clean and 
significant test of the behavior of soft radiation in QCD, and of the theoretical models 
employed in recent years to study power suppressed effects. 



2 The class of angularities 

We consider an e + e~ annihilation event with center-of-mass energy Q, generating a final 
state N, which we will take for now to consist of massless particles. The angularity, with 
weight a, of the state N, is defined as 

Ta(N) = -^J2Pi± e ~ ]ViKl ~ a) =^ y Eu i (&yO i ) a (l-\**O i \) 1 - a , (l) 

where Pi± is the transverse momentum of particle i relative to the thrust axis, rji is 
the corresponding pseudorapidity, r\i = lncot (0j/2), with 6i the angle with respect to 
the thrust axis, and Ui is the energy of particle i. The two definitions in Eq. (0) are 
equivalent for massless particles. 

The angularity distribution is defined as usual by summing over all final states, each 
weighed by its probability, according to 

^ = ^E \M(N)\ 2 ^r a -r a (N)), (2) 

where M(N) is the amplitude for the production of final state N. 

The thrust axis can be defined as the axis with respect to which Eq. $IJ is minimized 
at a = 0. The parameter a is adjustable in the range — oo < a < 2, with the upper limit 
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set by infrared safety. Angularity with a = is essentially 1 — T, with T the thrust pQ, 
while angularity with a = 1 corresponds to the jet broadening [3]. As discussed in |2T)] . 
for a > 1 recoil effects become important, so that the resummation of Sudakov logarithms 
must be modified, and consequently the pattern of power corrections changes, as pointed 
out for the broadening (a = 1) in jjj. In this paper we will be mostly concerned with 
the case a < 0, where power corrections are under good control. 

In the two-jet limit, r a — > 0, the distribution in Eq. (j2J) has large perturbative correc- 
tions containing powers of In (r a ), which have been resummed to all logarithmic orders, 
at leading power and for a < 1, in (201 • As usual, the resummation is performed at the 
level of the Laplace transform of the distribution, 



a ( v, a) 



dr n e" 



da 
dr a 



(3) 



At next-to-leading logarithmic (NLL) level, the resummation displays a nontrivial de- 
pendence on a, 



1 , r u 



1 I f du 

cr (z/, a) = exp < 2 / — 

0tot 1 J u 



v?Q 2 

1 



P\ 



A(a s (p ± ))(e- 1 - a ^r -l) 



^(a s (v^Q)) («• 



2/(2-a) 



(4) 



where and B(a s ) are the well-known anomalous dimensions acting as kernels of 

Sudakov exponentiation. The intricate a dependence of Eq. (HJ) simplifies at the level of 
leading logarithms, where one can easily invert the Laplace transform to find 



1 da 

°"tot dr a 



1 -a/2 vr 



a s ln(r a ) 
C F — - — exp 



1 



a/2 vr 



C F In 2 (r a 



(5) 



which displays a simple scaling with 1— a/2. This scaling of the perturbatively resummed 
cross section is however only approximate, and breaks down at NLL level, as can be seen 
in the explicit expressions given in the Appendix of Ref . [231 • 

Remarkably, the approximate scaling of the perturbative contribution with 1 — a/2 
is replaced at the level of leading power corrections by an exact scaling with 1 — a. 
To see this, one notes that the perturbative expression for the cross section given in 
Eq. (J3J) is ambiguous, due to the fact that the scale of the running coupling can vanish. 
As a consequence, at values of r a ~ A QCD /Q nonperturbative corrections must become 
dominant, and the perturbative expression needs to be supplemented by nonperturbative 
input to give a well-defined result. The structure of this nonperturbative correction 
can be deduced, following JSJ 123] > by introducing an infrared factorization scale k to 
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cut off the transverse momentum integration in Eq. (@J. The leading nonperturbative 
contribution arising from small transverse momenta can be evaluated by performing the 
integral over the Laplace variable u, keeping only terms scaling as powers of v/Q, while 
discarding terms suppressed by higher powers of the hard scale. The result can finally 
be written as a convolution of a perturbative contribution and a nonperturbative shape 
function, which in moment space is just a product, 



v 



In 



/a,NP 



OPT {y, K, a) / a , N p — , K 
i oo / \n 



1 + 



Q 2 - a ' Q 2 



n=l 



Q, 



(6) 
(7) 



The shape function / aj NP sums all power corrections of the form (v/Q) n , with unknown 
nonperturbative coefficients A n (fc), which can be formally expressed in terms of infrared 
moments of the cusp anomalous dimension A(a s ) as 



nn\ Jo 



dpi 



(plY /2 A(aM 



(8) 



As explicitly indicated in Eq. ©, terms of order u/Q 2 ~ a and of order /Q 2 have been 
neglected. At this level of accuracy, one finds that the only a-dependence of the shape 
function is through an overall factor 1/(1 — a) , which leads to the scaling rule [^H 1221 



0,NP 




(9) 



The derivation of the scaling rule in Eq. Q relies on two main assumptions. First, 
contributions from correlations between hemispheres are neglected, because the starting 
point is the NLL resummed cross section, which describes logarithmic corrections due to 
independent radiation off two back-to-back jets. In the more general resummed formula 
valid to all logarithmic orders [20], such correlations are present, but they contribute only 
starting at NNLL order. In addition, numerical studies indicate that inter-hemisphere 
correlations do not play an important role fOJEi]. One can furthermore argue that cor- 
relations between hemispheres due to particles whose decay products enter both hemi- 
spheres become non-negligible in the same range of the parameter a where also recoil 
effects need to be taken into account. The neglect of inter-hemisphere correlations is 
thus consistent with the resummation. The second assumption entering the derivation 
of Eq. (JUJ) is that nonperturbative soft radiation has the same properties under boosts 
as the relatively harder perturbative component. A success of experimental tests of the 
scaling rule would thus show that boost invariant dynamics dominates the differential 
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distributions at all scales and that coherent interjet radiation is nondominant in the 
relevant range of the parameter a. 

Of course, even if the above assumptions hold, there are further corrections present, 
suppressed relative to the dominant ones, as indicated in Eq. (JBJ). These corrections, 
which for thrust behave like v/Q 2 -, become important only in the extreme nonperturba- 
tive region 1 — T ~ (A QCD /Q) 2 , and are unlikely to play a role for phenomenology From 
a theoretical point of view, it may however be of some interest to compare the predictions 
of different models also for these subleading corrections. Considering the resummation, 
for example, one can keep terms neglected in Eq. (J0J), and parameterize them in terms of 
different integrals of the anomalous dimensions A(a s ) and B(a s ) in the infrared region. 
One finds that the pattern of subleading corrections can be characterized in terms of a 
subleading shape function as follows, 



a (u, a) = opt 0, k, a) / Q)NP ( ^, k) g a ^ P I — V 



(10) 



where 
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V 



\ Q 2 ~ a , 
Q 2 ) 



(11) 



The new nonperturbative parameters defining the subleading shape function g are given 
by expressions similar to Eq. (jHJ), 



A^(«,o) 



dp] 



, ., (pi)' 2 -"'" 2 A(a.(A)) 
nn\ Jo pj_ v ' v J 

^ 2 dp 2 



■r-B / , 1 f* dp Z ± / 2 xn/2 / 

nl Jo Pj_ 



(12) 



We notice that for both kinds of subleading contributions there is no simple scaling 
behavior with a. Furthermore, both contributions are suppressed at large v (that is, small 
T a ), with an increasing degree of suppression as a grows large and negative: specifically, 
subleading power corrections appear as functions only of the combination u/Q 2 ~ a , a 
feature that will also be found in the DGE formalism. 

In the following, we will study leading and subleading power corrections by means 
of DGE, and we will compare the results obtained with those arising directly from the 
resummation described above. 
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3 Dressed gluon exponentiation for angularities 



Dressed gluon exponentiation begins with a conventional renormalon analysis of the given 
event shape. One computes the characteristic function of the dispersive approach |12j . 
that is, the contribution to the cross section of a single gluon dressed with an arbitrary 
number of quark bubbles. This is referred to as the "single dressed gluon" (SDG) cross 
section. One proceeds by identifying the terms in the characteristic function which con- 
tribute to logarithmic behavior in the two-jet limit. These terms can be exponentiated, 
resumming all contributions where any number of dressed gluons is emitted without 
interfering. DGE reproduces resummation at NLL, provided the running coupling is 
defined in the bremsstrahlung scheme. Furthermore, all subleading logs are accounted 
for in the large /3q (large n/) limit. This way, DGE can detect the factorial growth 
of the coefficients of subleading logarithms, and provide methods to gauge the range 
of applicability of conventional Sudakov resummations. Furthermore, given an explicit 
representation of the singular behavior of perturbation theory, one may pick a definite 
prescription to deal with the resulting ambiguities. 



3.1 Angularities with a single dressed gluon 

Summing up infinitely many bubbles in the gluon propagator in the inclusive approxi- 
mation is equivalent to performing a calculation with an off-shell gluon with virtuality 
k 2 , and replacing the coupling with an effective, "timelike coupling" [30J. One then 
expresses the SDG cross section by 



1 da 



Ctot dr a 



^f-A^eQ 2 )^). (13) 
2pn Jo e v ' 



Here T is the derivative with respect to the gluon virtuality e = k 2 /Q 2 of the "charac- 
teristic function" J 7 , 

d 

F[T a , e) = -e— T{j a , e) , 

F(r a ,e) = J dx 1 dx 2 \M(x 1 ,x 2 ,e)\ 2 5(T a -T a (x h x 2l e)) . (14) 

The squared matrix element |M| 2 for the emission of a gluon of virtuality k 2 in the 
process 7* — > qqg, without coupling and color prefactors, is given by 

\w( \ |2 (xi + e) 2 + (x 2 + e) 2 e e 

x\ and x 2 , as usual, are the energy fractions of the quark and antiquark in the center- 
of-mass frame, and one can define X3 = 2 — x\ — x 2 , the gluon energy fraction. When 



7 



the gluon is "massive" , the limits of phase space are given by 



Xi + x 2 > 1 - e 
(l-x 1 )(l-x 2 ) > e. (16) 

The final ingredient in Eq. (jTSj) is the "timelike coupling" A eS , which is typically expressed 
in terms of a Borel representation as 




sin mi - 

A B (u) . (17) 



We now need to specify a suitable generalization of the massless definition of the 
event shape, Eq. (jTfl. for the case in which the emitted gluon has nonvanishing virtuality. 
Several observations are helpful in deciding how to perform this generalization. First, 
all logarithmic contributions to the SDG cross sections stem from the region where 
the gluon is either soft or collinear. Contributions from the region where the gluon 
is relatively hard do not give logarithmic enhancements, and the exact location of the 
boundary between the region of phase space where one of the quarks is dominant and the 
region where the gluon dominates is unimportant. We can thus concentrate on the phase 
space region where the thrust axis is the quark momentum, and X\ the largest energy 
fraction. Contributions in which the antiquark momentum dominates can be obtained 
by symmetry. Recall the expression for r a in the case of massless partons, in the region 
where the quark has the largest energy, 

T a (x U X 2 ) = ^ Ul - X 2 f- a l\x 1 + X 2 - If' 2 + (*! <-> X 2 )] . (18) 

X\ 1 J 

We would like our definition to reduce to Eq. (fTSj) as e — > 0. Furthermore, we require 
that for a = the definition should reduce to the massive definition of the thrust, as 
used for example in [THJ HZ1 UH1 EH] • In the phase space region at hand, this is simply 
t (xi,x 2 , e) = 1 — x\. Finally, as we will see, working analytically for generic a and with 
a massive gluon generates rather intricate expressions, so we must keep the definition as 
simple as possible in order to be able to perform the necessary integrations. 

Keeping these criteria in mind, we define the angularity with an off-shell gluon, in 
the region where the thrust axis is given by the quark momentum, as 

r.(xi,x 2) e) = (1 " Xl)1 ~ a/2 [(1 - x 2 - ef-^ix, + x 2 - 1 + ef' 2 

Xi L 

+ (a; 1 + x 2 -l + e) 1 - a / 2 (l-a; 2 -e) a / 2 ] . (19) 

Of course, other choices satisfying our criteria are possible. It can be shown, however, 
that different treatments of the gluon mass alter the value of r a by terms that are 
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suppressed by higher powers of the weight, roughly by a factor of T*~ a . Thus, they do not 
change logarithmically enhanced contributions. Once a definite prescription to include 
massive partons is chosen, the predictions within the DGE formalism are unambiguous. 
It should be kept in mind, in any case, that a comparison to experiment requires a 
detailed analysis along the lines of Ref. [2H] for the inclusion of hadron mass effects. 
To proceed, it is useful to change integration variables from xx,x 2 to 



c 



1 - X\ , 



1 + e 



1 — x 2 — e 



(20) 



In terms of these variables, the characteristic function can be written as 



1-2C 



+ 



(i-c + e ) 2 (i-o 

_C(i-C + <i + 0)(i + 
(£ + 2 (i-0 



e(l-C) 



C(i - C + c(i + 0)(i + 3 (i - C + c(i + 0) s 



e(l-C) 

C 2 (i + s 



5 r a 



c 



l-a/2 



i + e 



a/2, 



l-a> 



(21) 



Two comments are in order. First, the limits of integration do not correspond to the 
full phase space, but only to the region which generates logarithmically enhanced con- 
tribution, that is, the region in which the gluon is either soft or collinear to the quark. 
With a massive gluon, the collinear limit corresponds to x\ = 1 — e, x 2 = 0, whereas 
the soft limit is given by X\ = x 2 = 1 — y/e. These two values set the limits of the £ 
integration. The soft boundary of phase space for intermediate values of ( is given by 
x 2 < 1 — e/(l — xx), corresponding to the upper limit of the £ integration. The Sudakov 
region thus corresponds to values of £ close to the upper boundary of integration. The 
second observation is that, with this choice for angularity in the presence of a massive 
gluon, and in these variables, the 5-function defining the curves of constant angularity 
in Eq. (|2"Tj) does not depend on e, a feature that will be exploited below. 



3.2 Exponentiation and Borel representation 

The SDG cross section summarizes the probability that one dressed gluon is emitted. 
Upon summing over probabilities that many such dressed gluons are emitted indepen- 
dently, one obtains the exponentiated expression [13 IT8] 



-o (v, a. 



Ctot 



exp 



DGE 



1 da 



Ctot dr a 



(22) 
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Here we have extended the integration region beyond the support of the SDG cross 
section. This, however, does not change the result in the region of interest, at small r a 
or equivalently at large v. Inserting Eq. (|T3*|) and Eq. (fT7j) into Eq. (J2*2*|l . one may write 
down an explicit Borel representation for the exponent, although at this point one is still 
dealing with a five-fold integral. One finds 



In 



-a za a) 



Ctot 



9l 

2/3o 



du 



Q 2 



A 2 

v QCD , 



Sill 71U - , . 

A B {u) 

7TU 



J J dr a (e-^ - l) F{r a , e) 



o o 

(23) 

To exchange the order of integration, as done in Eq. ()23)1 . it has been necessary to 
evaluate the limiting values of angularity in the phase space region of interest, for a fixed 
gluon virtuality e. One finds that for any value of a the collinear limit of angularity 
is given by t°° = e + 0(e 2 ), while the soft limit is = sfe + 0(e) (notice that the 
corrections neglected here vanish for a = 0). As shown below, however, the exact form 
of these limits is not important: only the values of the energy fractions of the quark and 
of the gluon at the soft and collinear boundaries of phase space are relevant, and they 
are independent of the weight r a , and solely dependent on the kinematics. 

The advantage of writing the exponent as in Eq. (}2"3*|) is that it is now possible, before 
performing phase space integrals, to take the derivative with respect to e of Eq. (|2*Tj) . 
thanks to the fact that the 5-function defining r a does not depend on e with our choice 
of variables. One may then use the 5-function to perform trivially the r a integral. Dis- 
carding terms that do not contribute Sudakov logarithms, we can write the Borel repre- 
sentation of the exponent as 



In 



1 



Ctot 



■cr(z/, a) 



Cf_ r du ( Q 2 ' 

2/Wo AAqcd, 



with the Borel function B(u, u, a) given by 

oo V~e ( CA-1 ! 

B{u,v,a) = Jdtt- u Jdc\ J 



C 2 (i + 2 



sin iTU 



A B {u) B(u, v, a) 



( fa/2 

exp -uC, {2 - a)/2 5 




l-a\ 



(24) 



.(25) 



The factors of —1 in Eq. ()25j) are due to virtual corrections. They contribute terms 
independent of a, and are thus identical to the thrust EH]- They do not need to 
be considered anew. We will now concentrate on the contributions associated with real 
gluon emission, labelled in the following by an additional subscript, Br(u, u, a). 
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It is possible to arrive at a one- dimensional integral representation of B R (u,u,a). 
First we perform the integration over £ in the first term of Eq. ()25|). and change variables 
from ( to £ = C/e— 1 in the second term. Finally, integrating over e results in the following 
expression 



B R (u, is, a) 



oo 



2 1 

+ 



1 



l + £ (i + 2 i-«(i + O s 



2m 



2 -a' 



2m 
2 - a' 



a/2 



l-a 



2-0/2 



(26) 



where 1.F1 is the confluent hypergeometric function, also known as Kummer's function 
of the first kind, defined by [3T] 



t F! (a; /?;*) = £ 



fro 09)* W ' 



(27) 



where the Pochhammer symbol is defined as 



[a)k = a(a 4- 1) . . . (a + k — 1) 



T(a + k) 
T(a) 



(28) 



Clearly, Eq. ()26|) cannot be evaluated directly. Nevertheless, in the physically interesting 
limits, the soft and collinear region of phase space, results in closed form can be obtained, 
and from these results the form of power corrections can be inferred. Before discussing the 
results for general a, however, we pause to note that we can test Eq. ([SB} by considering 
the case of the thrust, a = 0, where the answer is known [T7| I18j. For a = 0, indeed, 
Eq. ()26|) can be explicitly integrated. From the resulting expressions one can see that 
only the factor 2/(1 + £) in the square bracket of Eq. (|2*Bj) contributes in the soft limit, 
whereas the remaining contributions are of purely collinear origin. The final result for 
B R (u, v, 0) is 



B R (u, v, 0) = v 2u T(-2n)- - v u T{-u ' " 

u 



1 



1 



u 



u 



u 



(29) 



up to corrections of order 1/v. Our result is in complete agreement with Refs. fT7\ I18j. 
Soft contributions produce singularities at half-integer values of u, corresponding to 
power corrections of the form (A QCD u / Q) p '; collinear contributions have poles at integer 
u, giving power corrections behaving as (Aq Cd z//Q 2 ) p . We note also that the result in 
Eq. (|2~9"j) is infrared safe at leading power. In other words, the poles at u = (corre- 
sponding to Q°) cancel between the soft and the collinear contributions. We will verify 
this key property of infrared safety also for general a. 
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4 The pattern of exponentiated power corrections 



4.1 Power corrections from soft radiation 



We now turn to the soft contribution for general a, which is the source of the scaling in 
Eq. (jUJ). In the region of phase space corresponding to soft emission one has ( ~ y/e, 
so that £ is large, behaving as In this region, we will first evaluate Eq. (|2T)j) by 

making use of an asymptotic expansion for the Kummer function. Then we will go back 
to the characteristic function T{r a , e) and evaluate it in the same limit, to trace the 
contributions that display scaling back to the matrix element. 

To treat the soft limit, we can use the following integral representation of the confluent 
Kummer function [3*2] . 



iFAa; (3, z) 



a) Jo 



I 1 dte zt t a - l (l-t) 
Jo 



-a+P-l 



(30) 



We then see that the complicated argument of the Kummer function in Eq. (jzbj) appears 
only in the exponent, and can be expanded for large £ and u, with z//£ held fixed. One 
finds 

i+e- a 



cxp 



-v 



t ^a/2_ 





ut 


= exp 





+ 



(i + 2 ~ a/2 . 

The soft contribution then simplifies considerably, and becomes 



(31) 



^° ft («, v, a) 



CXJ 1 



x 



o 
2 



+ 



i + £ (i + o 2 i-u(i + o 



+ o[- 

.V 



(32) 



Since we are interested in the large £ limit, clearly at leading power only the first term 
in the bracket of Eq. (}3*2*|) contributes, while the remaining two terms are suppressed by 
powers of \jv. The integral can be performed, and one finds that the final answer for 
the soft contribution is 



1 



,2m 



-2u) - + 
u 



(33) 



Comparing this with the case for a = 0, Eq. (|29|h we see that we recover exactly the 
prediction of Eq. Q. We have thus confirmed the result of j^S], obtained there by 
analyzing the resummed expression. 
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4.2 An alternative derivation of the scaling rule 



To perform the integrations leading to Eq. (}33|) for generic a, it has been necessary to 
make use of some special function technology. This technology may obscure the physics 
underlying the simple scaling, and indeed the result of Eq. (|33|) may appear somewhat 
surprising, in view of the fact that, for example, Eq. ()32j) still has an overall factor 
of 2/(2 — a), corresponding to the scaling of leading logarithms but not to the scaling 
of power corrections. It is thus instructive to recover the result by going back to the 
characteristic function, Eq. (|21j) . and identifying the origin of the terms responsible for 
the leading soft contribution. 

The first observation to simplify Eq. ()21|) concerns phase space. As observed before, 
the Sudakov region corresponds to a neighborhood of the upper limit of the £ integration 
range, whereas the lower limit does not contribute to logarithmic enhancements and can 
be changed. To simplify the calculation, we can for example take as a lower limit £ = 0, 
as was already done in Eq. (|25|). The expression for T becomes then slightly more 
manageable by changing variables from £ to u = 1 + £. One finds 



du J d( 

1 6LU 



i-C 



;i-c + e ) 2 + (c + ^-i) 2 



c^ 2 



w(l-C + ew) u 3 (l-( + euj) 
5 (r a - C W2 /aH) , 



where we define 



a (oj) = (uj - 1) 



l-a 



(34) 



(35) 



Notice that fo(u) = 1, which of course simplifies things considerably for the thrust. For 
any a, the ( dependence in the 5-function is sufficiently simple to be used to perform 
the ( integration, yielding a rather cumbersome expression, which however has some 
interesting features. One finds 



F(T a , e) 



W , 
2 1 7 du 



2 Ob Tn J UJ 

1 



1-Ca 



where 



Ca 



( a u (1 - C. + euf 



2 

2-a 



_(l-( a + eu) u\l-C a + eu) 

e(c a -eu) e(V~e-c a ) , (36) 



(37) 



\faHj [1 + (CO - l)l-«] — 

Several observations help disentangling Eq. (|36|). First, the second ^-function is not 
relevant to the logarithmic behavior and can be neglected. Its effect would be to split 
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the integration range for oj in two subintervals, however all leading contributions come 
from neighborhoods of the upper and lower limits of integration. A second important 
point is the fact that the leading singularity in r a is now explicitly factored out. One can 
then evaluate the integral to leading power in r a . The constraint imposed by the first 
^-function can also be considerably simplified. The integral, in fact, has support on the 
region defined by 

Tn . ( OJ W 2 



7^^{— ) i 1 + J • (38) 

This constraint cannot be solved exactly for uj. One notices, however, that in the physical 
region, e < r a < y/e, and for a < 0, the left hand side of Eq. (|3*5j) is a parametrically 
large number. Within the integration range, on the other hand, the right hand side 
becomes large only near the boundaries, as uj — ► oo or as uj — > 1. One can then solve the 
^-function constraint in these two limits, obtaining respectively 



and 



uj < UJ + 



UJ > UJ- 



-l-a/2 



A-2/a 



1 + 



Ta 



■2/a 



OJ 



UJ 



OO 



1 . 



(39) 



(40) 



The result of these manipulations is that the integration region in Eq. (jBEj) shrinks at 
both boundaries, with both integration limits now dependent on r a . Since soft contribu- 
tions arise from the region of large uj, we will concentrate on Eq. (|39jl. although similar 
arguments could be used with Eq. (|40|l . This will be sufficient to recover the scaling rule. 

To complete the calculation, we must now approximate the matrix element. To this 
end, note that for large uj one can approximate Eq. (J37J) by 



(41) 



One sees that ( a is a small number in the relevant range of a, and we can expand 
the integrand of Eq. (|3f)jl around ( a = 0. The resulting soft approximation of the 
characteristic function is 



F(r a , e) 



I soft 



2 - a T n 



UJ 



1 



OJ 



(42) 



The integration can now be performed and compared with the result for the thrust. 
At large uj the dominant contribution clearly comes from the first term in the square 
bracket. Comparing with the case a = 0, we see in fact that, while all terms contribute 
to the Sudakov limit, the last three terms give sub leading corrections associated with 
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collinear radiation. This remains true for generic a, as noted in the following subsection. 
Concentrating on the first term, we can finally see how the scaling rule arises in the 
context of DGE. The only term in the cross section contributing in the soft limit is, as 
might be expected, the logarithmic integral over the bremsstrahlung gluon spectrum. 
The gluon mass then acts as an infrared cutoff on this integral. The power of e in the 
upper limit of integration u + is precisely the one required to cancel the overall factor of 
2/(2 — a), and to replace it with the scaling factor of Eq. Q. This precise power arises 
uniquely from the definition of angularity, as expressed by the 5-function in Eq. (}3"4"j) . an d 
can easily be traced back to the exponential weight given to pseudorapidity in Eq. ([TJ. 
The final result in the soft limit is simply 

nr a ,e)\ soft = --^-^He), (43) 

1 CI T a 

where sub leading collinear contributions and terms independent of e, which do not con- 
tribute to the logarithmic behavior, have been omitted. Power corrections of collinear 
origin can be treated similarly, giving results consistent with the ones outlined in the 
following subsection. 



4.3 Power corrections from collinear radiation 

At this point we have various methods at our disposal to analyze power corrections of 
collinear origin. We can for example complete the analysis of Section 14.21 including 
the effects of the lower limit of integration, and then inserting the results in the Borel 
exponent. Alternatively, we can go back to Eq. (J2T?|) and study it in this limit. Collinear 
power corrections arise at large u, but small £. In this limit, Eq. (j2Sj) yields 



2-o V 2-o 



oo 

2u 





"2 2 1 
x - — + 



l+£ (l + O 2 i-«(i+O s 



+ <?(-)■ (44) 



v 



We see immediately that collinear contributions are suppressed by non-integer powers 
compared to the leading soft piece, as found above in the derivation of Eq. 0. This 
non- integer power behavior can also be seen from Eq. (|44jl . which has singularities at 
non-integer u due to the T-function, and further singularities due to the integration over 
£, which can be readily performed for any specific value of a. For generic a, the leading 
poles on the positive real u axis can also be determined by studying the integrand of 
Eq. (J44j) in the limits of large and small £. Determining the full analytic structure of 
Eq. (J44j) as a function of a is however considerably more difficult. In the case when a is a 
rational number, we have been able to obtain a closed expression in terms of generalized 
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hypergeometric functions p F q . The detailed expressions are listed in Appendix A. Since 
all methods give consistent results, we will just summarize here the structure of poles 
on the positive real axis of the Borel variable u, and we will outline the corresponding 
pattern of power corrections. 

First, let us note that our complete result for the Borel function B is infrared safe, 
as it must. In fact, as was the case for the thrust, the poles at u = cancel between soft 
and collinear contributions. To see this, note for example that the intricate expression 
in Eq. (J58J) simplifies at u = to 

sr(u ~ „,„,„) = ^«r(-^- M ) ^-r (-?fi^») , (45) 

where we have suppressed terms nonsingular at u — 0. This clearly cancels against the 
contribution of Eq. at u — 0, as expected. 

Next, we note that Eq. (J44]) appears to have an explicit pole at u — 1, which would 
correspond to a correction of order O (u(A QCD /Q) 2 ). This pole, however, is cancelled 
by the explicit factor of shi7ru in Eq. (J23j) . as was the case for the thrust (where power 
corrections of this form are present only as an effect of the T-functions in Eq. (J29j) ) . 

The general structure of poles in u, for rational a, can be deduced from Eq. (J55j) in 
Appendix A. There are several infinite towers of poles. From Eq. f)44j) we can directly 
read off the first pole in each tower. We find singularities at 

2 - a 



u = u 2 = , (46) 

a 

12-a 
2 1 — a 

Recall that a < in the range of validity of our approach. The crucial fact here is 
that all these singularities are accompanied by a common factor of u 2u ^ 2 ~ a \ as seen in 
Eq. (|44jl . When combined with the locations of the poles in Eq. (|46p. we see that all 
collinear power corrections are expressed in terms of a single combination of v and of 
the scale Q, namely h , (A QCD /Q) 2 ~ a , possibly further raised to a non-integer power. This 
result agrees with the estimate extracted for the resummation, Eqs. (|l() jl -pi jl , although 
the detailed pattern of subleading singularities is different. We conclude that collinear 
power corrections are suppressed as predicted from resummation, and are expected to 
become important only for extreme values of the angularity r a , r a ~ C(A QCD /Q) 2_a . 
The suppression grows as a becomes large and negative, although numerically the effect 
could be partly compensated by the fact some of the corrections may be further raised 
to small non- integer powers, of order 1/(1 — a). 
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It is also important to note that the intricate structure of towers of subleading poles 
arising from Eq. (J55|) is not as model-independent as the leading poles connected to 
soft radiation. For example, a different choice for the massive definition of angularity, 
Eq. (JH1J), might kinematically generate non-integer power corrections of comparable size. 
We emphasize, in any case, that the leading structure expressed by Eq. is unaffected. 
Furthermore, we believe that the parametric dependence of collinear power corrections 
on the ratio u(A QCU /Q) 2 ~ a , which is found via resummation, and confirmed by our DGE 
analysis, is a stable feature, suggesting that the scaling rule should hold with increasing 
precision for negative a. 

5 Conclusions 

In this work we have verified that the universal scaling of the leading power corrections 
within the class of angularities, discovered in Ref. |2Sj from soft gluon resummation, 
is unchanged when, in addition to large Sudakov logarithms, also renormalon enhance- 
ments are taken into account. Furthermore, we have been able to determine the form of 
subleading corrections within a specific scheme to treat massive partons. These correc- 
tions, of collinear origin, are suppressed by non-integer powers of the hard scale relative 
to the leading ones that originate from soft radiation. Different models of these collinear 
corrections further suggest that the relative suppression grows as the parameter a be- 
comes large and negative. A detailed phenomenological study along the lines of Refs. 
[T7| [THj could help to further constrain the nonleading effects, and would be necessary 
for a comparison with experiment. We defer such a study to future work, due to the lack 
of corresponding experimental data. 

Assuming subleading corrections are negligible, as indeed our results imply, the scal- 
ing allows to predict the distributions of the whole class of angularities, including non- 
perturbative corrections, in terms of a single shape function, which can in principle be 
determined from data by considering a specific choice of the parameter a, for example 
the thrust, a = 0. Since there are no free parameters, an experimental determination 
of distributions for other values of a would certainly give valuable information about 
the properties of nonperturbative corrections. Also, since the perturbative cross section 
has a scaling behavior with a different from the non-perturbative shape function, com- 
parison of theory and experiment might help to disentangle corrections due to missing 
higher-order perturbative information from power correction effects. An experimental 
test of scaling would determine to what extent the boost invariance of soft radiation in 
the two-jet limit, which is established at the perturbative level, also applies to nonper- 
turbative effects. We hope therefore that an analysis of experimental data for the class 
of angularities will be made in the near future. 
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A Evaluation of the collinear contribution 



Here we sketch the evaluation of the collinear contribution to the DGE cross section for 
rational a, that is, we consider the case 

P 

1 — a = -, p,q positive integer, p > q. (47) 
Q 

This is sufficient for our purposes, since any irrational number can always be well 
approximated by a rational one. 

The collinear limit is equivalent to the limit in which v —>■ oo faster than £ (although 
£ varies between and oo). In this limit the collinear contribution is given by Eq. ([44)1 . 
In the following we will abbreviate Eq. (J4*4*|) by 

B^iu, v, a) = -l-r (--^-) v& £ c m J™, (48) 
/, a \ z a j m=1 

where 

oo 

r = y c /a^"(i + 0~^ ,1 " m (i + ^ 1 " a )^ r (49) 



d = 2, c 2 = -2, C3 = _!_. (50) 

1 — U 

We use the Mellin-Barnes representation of (1 + £ 1 - a ^ 2u /@- a ) j a — \ — pjq^ 
v s ; 2m Jc r (—2*0 



The contour C runs along the imaginary axis, to the left of Re(ct) = 0. With this we 
can rewrite Eq. (j4T?j) in terms of p and q as 

^ = r da r(-a)r(a-^ M )r(m-i-ga)r(i + ^ + ;q) 

We change variables from a to 5 = a/q, with C the contour in the new variable, and 
use the following properties of the T-function [21 E2] : 

T{nz) = n nz - 1/ \2nf- n]/2 l[r(z + -) , n integer, (53) 



fc=0 



Res Q {r(a- &)/(<*)} = £^-/-/(a = 6-n), (55) 
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where in (|55|) the residues of the r-function are taken with respect to a. This allows us 
to rewrite as 



x nrfa-— w + nrfa + 7 ^-« + — V (56) 

We now close the contour to pick up the residues at a = — ~ q — n, n integer, where 
j runs from to q — 1, and a = — — j^§w, u ~ n i where h runs from to p — 1. We 
thus have p + q residues, evaluated according to (J55j) . After a bit of algebra, using Eqs. 
(153(1 and (|54j). as well as the definition of the Pochhammer symbol (J28j) . we can rewrite 
the result in terms of generalized hypergeometric functions, which are defined as pH] 

p 

p F q ({a u ... } a p };{f3 1} ... } (3 q };z) = J2 — • ( 57 ) 

B =° n! n (&)» 
3=1 



In order to find the poles in tt with respect to a, it is instructive to partially rewrite the 
result in terms of a where possible (recall that p — (1 — a) q): 



r=0 i=0,jVr 



Q 



x (2^) 9 F( 2 -a) 9 -i(M};{^};(-l) (2 - a)(? 



2-a' 1 



r (-^) 



ri 


^(1 — a) r + m — 1 — 2 ^_ °^ ttj 


ir(l + «-(l-a)r) 


r (m + ^ 


l 



(1— a)q— 1 (1— a) g— 1 / • , \ 

+ £ ^a-.)-/^)^ p r to- (58) 



X ( 2 -a) g F(2- o)9 -l (Kl; {B'}; (-l)^ 



r(™ + ^) r (-^«) 
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with the following arguments of the generalized hypergeometric functions 
{A} = [a 1; . . . ,a q ,6>i, . . . ,a ( i_ a)g } 

{B} = {p 1 ,-..,P^ r+1 ,...,P q Ji,...J(i- a)q } (59) 

2 r - 1 i 

(2-a)g g g 

2 r m — 2 i 

(2 — a)g g (1 — a) g (1 — a)g' 

r + 1 i 

A = 1 + i^r + l, 

q q 

~ „ u r i 

A = l-7i ^ + 



(1 — a) g g (1 — a) g' 
{.A'} = {a' 1 ,...,Q!j,5' 1 ,...,5( 1 _ a)(f ] 

{#} = {P[,...,^J[,---Ji\ ¥t+1 ,---Jlx-a) q } (60) 



, a u t + 1 1 i 

a- = 1 1 — , 

2 — a (1 — a) g (1 — a) g g g' 

a u t + m — 1 z 



a = z ^ + ^ ^ + 



2 — a (1 — a) g (1 — a) g (1 — a) g' 

„ t + 1 it 1 z 

A = 1 + 7^ r- + 7; r- + , 

(l-a)g (l-a)g g g 

ft = 1+n r »^ + L 

(1 — a) g (1 — a) g 

Eq. ()58|) inserted into Eq. (jlHj) gives the final answer for the collinear contribution. For 
a = or equivalently p — q — 1, we reproduce the collinear part of of course. 

The poles in w can be read off from (J58|) . using the properties of the Gamma- and 
hypergeometric functions involved (211 E2]- The T-function T(z) has simple poles at 
z = —n, n integer, with residues (— l) n /n\, see Eq. (J3SJ), the hypergeometric p F q (J57|) 
has simple poles only in (3j at (3j = —n, n integer, with residues 

(-l) n - 

p F q ({ai} Pj-i, -n, p j+1 , ...,/3 q };z), (61) 



n! 

where p F q denotes the regularized hypergeometric function 

v 



oo n II \(%i)n 

z 



,F q ({ ai , . . . , a p } ; {A, • • • , A,} ; ^) = E -T • ( 62 ) 

«=o n - n r(n + A) 

20 



This means, that despite its appearance, (|58j) has poles in u only in the T-functions, and 

2p 

in the fit and Note that any pole in u at u = p is accompanied by a factor of z/ 2 ~ a , 
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